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We use appropriately defined short ranged reference models of liquid water to clarify the different 
roles local hydrogen bonding, van der Waals attractions, and long ranged electrostatic interactions 
play in the solvation and association of apolar solutes in water. While local hydrogen bonding in¬ 
teractions dominate hydrophobic effects involving small solutes, longer ranged electrostatic and dis¬ 
persion interactions are found to be increasingly important in the description of interfacial structure 
around large solutes. The hydrogen bond network sets the solute length scale at which a crossover in 
solvation behavior between these small and large length scale regimes is observed. Unbalanced long 
ranged forces acting on interfacial water molecules are also important in hydrophobic association, 
illustrated here by analysis of the association of model methane and buckminsterfullerene solutes. 


I. INTRODUCTION 

Hydrophobic interactions play a key role in phenomena 
ranging from biological processes like protein folding and 
membrane formation to the design of water-repellent ma¬ 
terials m ■ Thus, significant effort has been devoted to 
studying the behavior of apolar moieties in water. In pi¬ 
oneering work, Stillinger argued that hard sphere solutes 
smaller than a critical radius Rc can be inserted into liq¬ 
uid water while maintaining the hydrogen bond network, 
but for solutes with a radius larger than Rc bonds must 
be broken, generating a molecular scale interface with 
properties resembling that of the liquid-vapor interface 
in water [J]. More recent work has confirmed the basic 
features of this idea and put the arguments on a firmer 
statistical mechanical foundation MEIE]. 

While this qualitative description of the length scale 
dependence of hydrophobic hydration seems physically 
very reasonable, it focuses only on the hydrogen bond 
network of water and makes no mention of the van der 
Waals (VDW) attractions and long ranged multipolar in¬ 
teractions between water molecules or of the VDW at¬ 
tractions that would be present between a more realistic 
solute and the solvent. Moreover, a qualitatively sim¬ 
ilar length scale transition is seen in a dense Lennard- 
Jones (LJ) fluid near the triple point, with the forma¬ 
tion of a “dry” vapor-like interface around a large hard 
sphere solute [6j. In that case clearly there are no hy¬ 
drogen bonds and the transition is generated solely by 
unbalanced VDW attractive forces arising from solvent 
molecules far from the solute. 

Consideration of such unbalanced forces is an essential 
ingredient in the Lum-Chandler-Weeks (LCW) theory of 
hydrophobicity [5], which uses the same basic framework 
to describe hard sphere solvation in simple liquids and in 
water, differing only in the thermodynamic parameters 
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needed as input to the theory p 0 El • Indeed LCW the¬ 
ory has been criticized for not treating hydrogen bonds 
and other distinctive features of water more explicitly 
and there has also been considerable debate about pos¬ 
sible effects of solute-solvent LJ attractions on the pro¬ 
posed length scale transition in water [3j. Thus it seems 
useful to explore in more detail the varying roles hydro¬ 
gen bonds, VDW interactions, and long ranged multi¬ 
polar interactions play in hydrophobic solvation, and to 
determine what analogies exist to solvation in simple, 
non-associating fluids. 

To that end, we build on our previous work <|9] us¬ 
ing truncated water models pm, and exploit the under¬ 
lying ideas of perturbation mm and local molecular 
field [5; |T3] (LMF) theories of uniform and nonuniform 
fluids, respectively, to study hydrophobic solvation and 
association from small to large length scales. We em¬ 
ploy short ranged variants of the SPC/E water model to 
show that small scale solvation and association in water is 
governed by the energetics of the hydrogen bond network 
alone. However when the solute is large and the hydro¬ 
gen bond network is broken at the hydrophobic interface, 
water behaves in a manner qualitatively similar to a sim¬ 
ple fluid, with unbalanced LJ attractions dominating the 
solvation behavior. 

In the next section, the truncated water models are 
briefly introduced and our simulation methods are de¬ 
tailed. Section III examines the roles of unbalanced dis¬ 
persion and electrostatic forces in determining the equi¬ 
librium solvation structure around small and large apo¬ 
lar solutes. The strength of the hydrogen bond network 
around small solutes is then analyzed by perturbing the 
hydration shell in Section IV. The role of this network 
in setting the length-scale for the crossover in solvation 
thermodynamics is then studied in Section V. The origin 
of entropy convergence is briefly discussed in Section VI. 
Finally, the hydrophobic association of model methane 
and fullerene molecules is studied in Section VII. Our 
conclusions and a discussion of the implications of this 
work are given in Section VIII. 
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II. MODELS AND SIMULATION DETAILS 

Hydrogen bonds in most classical water models arise 
from “frustrated charge pairing”, where an effective pos¬ 
itive charge on a hydrogen site of one molecule tries to 
get close to a negatively charged acceptor site on a neigh¬ 
boring molecule 0. This strong attractive interaction is 
opposed by the overlap of the repulsive Lennard-Jones 
(LJ) cores and the presence of other hydrogen sites in 
the acceptor molecule. As a result, short ranged versions 
of the full water model where Coulomb interactions are 
truncated at distances larger than the hydrogen bond 
length and with only truncated LJ core interactions if 
desired can still give a very accurate description of the 
hydrogen bond network and pair correlation functions in 
bulk water mm- 

In this work, we use the extended simple point charge 
(SPC/E) model of water [Hj and two previously devel¬ 
oped short ranged variants of this model 0 to examine 
hydrophobic hydration and association as the solute per¬ 
turbs the hydrogen bond network. The truncated mod¬ 
els provide a hierarchical framework for disentangling in 
such classical models the separate contributions of (i) 
strong short ranged interactions leading to the hydrogen 
bond network, (ii) longer-ranged VDW attractions be¬ 
tween water molecules and with the solute, and (iii) long 
ranged dipolar interactions between water molecules. 

The first such model, the Gaussian-truncated (GT) 
water model, has full LJ interactions but truncated 
Coulomb interactions 0 m- It thus lacks the long 
ranged electrostatics necessary to provide a description 
of the physical multipolar interactions that act over 
large distances. We also utilize the Gaussian-truncated 
repulsive-core (GTRC) model 0, where both the long 
ranged electrostatic interactions and the long ranged LJ 
attractions have been removed. The GTRC model gen¬ 
erates a minimal reference network model that captures 
very well the structure of the local hydrogen bond net¬ 
work in bulk water while ignoring effects of the remaining 
long ranged Coulomb and dispersion interactions. 

In order to compare the SPC/E water model at a 
pressure of P = 1 atm with the short ranged GT and 
GTRC models in the work presented below, the latter 
two models were simulated at corrected pressures yielding 
the same density using the pressure corrections described 
earlier mm- There it was shown that simple analytical 
corrections to the pressure can bring the bulk densities of 
these three models into quantitative agreement. All data 
presented in this work were obtained from molecular dy¬ 
namics simulations performed in the isothermal-isobaric 
ensemble (constant NPT) using a modified version of the 
DL_POLY2.18 software package El- Constant temper¬ 
ature and pressure conditions were maintained through 
the use of a Berendsen thermostat and barostat, respec¬ 
tively |T8] . The evaluation of electrostatic interactions in 
simulations of the full SPC/E water model employed the 
Ewald summation method [Sl¬ 
it is instructive to compare the solvation behavior of 


water to that of a simple LJ fluid at an analogous state 
point throughout this work. Therefore, following the 
work of Huang and Chandler 0, we also study a LJ fluid 
at a state point near the triple point, where the potential 
is truncated and shifted at 2.5cr. This LJ fluid is studied 
at a reduced temperature and pressure of T* = ksT/e = 
0.85 and P* = Pa 3 /e = 0.022, respectively, correspond¬ 
ing to a reduced density of p* = pa 3 = 0.70. In order to 
study the analogous short ranged reference fluid, we use 
the same repulsive force truncation of the LJ potential 
as was used for the GTRC water model, and study the 
model at a mean-field corrected pressure that accounts 
for the lack of LJ attractions 0. 

We should emphasize that the above-mentioned short 
ranged GT and GTRC models are not being used in this 
paper as replacements for standard long-ranged models 
such as SPC/E or to give accurate representations of 
most properties of real water. Rather, we utilize these 
models as analysis tools to examine the different roles 
the hydrogen bond network as described by the GT or 
GTRC models, long ranged dispersions, and dipolar in¬ 
teractions play in determining the properties of systems 
containing liquid water. 

However, the GT model describes very well pair cor¬ 
relation functions and hydrogen bond statistics in bulk 
water, and as we discuss further below, it also captures 
many features of the water density in nonunifornr envi¬ 
ronments including the basic length scale transition for 
hydrophobic solutes 0. But thermodynamic and par¬ 
ticularly electrostatic properties depend sensitively on 
the long ranged Coulomb interactions and GT results 
need corrections for quantitative accuracy. Acharya and 
Garde 05j have recently carried out a detailed study of 
the strengths and weaknesses of the GT model as a sim¬ 
ple water model in a variety of settings, including both 
hydrophobic and ionic solvation. 


III. THE INFLUENCE OF LONG RANGED 
INTERACTIONS ON INTERFACIAL 
STRUCTURE 


In this section, we examine the role of the various un¬ 
balanced forces in determining the interfacial structure 
of water near a hydrophobic solute. The solute is con¬ 
sidered to be a uniform density of LJ particles, such that 
its interaction with water can be represented by an inte¬ 
gration of the LJ potential over the volume of the solute, 
resulting in the integrated “9 — 3” potential of Huang and 
Chandler J2T] 
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FIG. 1. Density distributions around solutes of radii Hhs ~ 
2 A (a) and Rus ~ 20 A (b). The inset in (b) depicts the 
renormalized portion of the LJ LMF for GTRC water in units 
of k B T. 


where r± = r ± Rs■ The parameters of the poten¬ 
tial are chosen to mimic paraffin, such that the den¬ 
sity of LJ sites, energy, and length scale are given by 
p = 0.0240 A -3 , e sw = 0.882 kJ/rnol, and a sw = 3.468 A, 
respectively [21] . Furthermore, in order to make this par¬ 
ticle as hydrophobic as possible, only the repulsive part 
of the potential is used, such that the solute-water in¬ 
teraction potential used in the MD simulations is given 
by 


Go,sw(^ ) 


U sw (r) - U sw (r 0 ), r < r 0 
0, r > r 0 


where ro is the location of the minimum of the poten¬ 
tial. Finally, we should note that the size of the particles 
is better represented through an effective hard-sphere ra¬ 
dius, i?HS> rather than the size parameter Rs found in the 
potential. This effective radius can be estimated as [22] 

/»oo 

Rhs ~ dr {1 - exp [-/3U 0 , SW (r)]} , (3) 

Jo 

where (3 = (k B T) , and will be reported as i?ns herein. 

The hydration structure around small solutes has been 
postulated to be a direct consequence of the need for wa¬ 
ter to maintain its hydrogen bond network. A small so¬ 
lute can be “inserted” into bulk water with the network 
continuing around the solute without breaking hydrogen 
bonds. Indeed, in the small solute regime we find that the 
nonuniform densities of GT and GTRC models around 
an apolar particle are nearly identical to that of the full 


SPC/E water, dramatically confirming that local hydro¬ 
gen bonding dictates the hydration structure in this limit 
(Figure [T}i). 

In the large solute limit, Figure[lJ>, the density profiles 
of SPC/E and GT water are still very similar, demon¬ 
strating that long ranged electrostatic interactions have 
an almost negligible influence on this measure of inter¬ 
facial structure. GTRC water, on the other hand, has a 
p(r) markedly different from that of SPC/E water. 

Removal of the LJ attractions from the bulk liquid in 
GTRC water eliminates the phenomena of drying, and it 
evidentially wets the surface of the solute. According to 
LMF theory m , we can account for the averaged effects 
of the neglected LJ forces by using a renormalized solute 
field 

</R J ( r ) = U 0iSW (r)+J dr' [p R (r') - p B \ m (|r - r'|), (4) 

where quantities obtained in the presence of the effec¬ 
tive field are indicated by the subscript ‘R’ throughout 
this work, p B is the bulk density of the fluid, and u\{r) 
is the attractive portion of the LJ potential. The use 
of this renormalized field recovers drying behavior and 
brings the density profile of GTRC water into qualita¬ 
tive agreement with that of the SPC/E and GT mod¬ 
els, as illustrated by the curve labeled ‘GTRC-LJLMF’ 
in Figure [Tj >. The renormalized portion of the LMF, 
</>p/(r) = <^ J (r) — L r o jSW (r), provides an effective force 
that pushes solvent molecules away from the solute, as 
shown in the inset of Figure [l]r. 

From the data presented in Figure [lj we can conclude 
that the unbalanced forces arising from LJ attractions 
are the driving force for drying at extended hydrophobic 
interfaces. Indeed, we have previously shown that the net 
force on a water molecule at an extended hydrophobic in¬ 
terface from long ranged electrostatics is much smaller 
than that from LJ attractions 9j. Nevertheless, long 
ranged electrostatics play a subtle but important role in 
determining the orientational preferences of water and 
properties dependent upon this orientational structure. 
One such quantity is the electrostatic or polarization po¬ 
tential <I)(r) felt by a test charge 

<F(r) = — f dr'£(r') 

Jo 

r r dr' r' 

= -/ -Ta/ drVVtr"), (5) 

Jo r Jo 

where p q (r) = (JT qiS(r — r*)) is the ensemble averaged 
charge density of the system and £(r) is the electric field 
due to the polarization of water molecules induced by the 
presence of the solute. 

The polarization potential of SPC/E water, shown in 
Figure Sb reaches a constant value of approximately 
500 mV in the bulk region, consistent with previous 
determinations of interface potentials at extended hy¬ 
drophobic interfaces for this water model m■ Removal 
of the long ranged electrostatic interactions in GT wa¬ 
ter leads to an approximate charge density that does not 
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predict this plateau in the bulk region, Figure [2jr. Thus 
there is a net electric field £(r) in this system, even far 
from the solute surface as shown in Figure [2]:). The ap¬ 
pearance of a non-vanishing electric field in the bulk of 
GT water is associated with an over-orientation of inter¬ 
facial OH bonds toward the solute surface. This is evi¬ 
denced by a larger peak at $oh ~ 0° in the probability 
distribution P ($oh) for interfacial GT water molecules 
in comparison to that observed for SPC/E water, shown 
in Figure [2};, where $oh is the angle formed by the OH 
bond vector and the oxygen-solute vector 

The increase in the number of OH groups pointing to¬ 
ward the interface in GT water is driven by the tendency 
to maintain the hydrogen bond network alone. This re¬ 
sults in the formation of an overly ordered dipole layer 
at the interface, demonstrated by the peak at 9 M ss 60° 
in P(9fj), shown in Figure where 9^ is the angle 
formed by the dipole vector of water and the oxygen- 
solute vector. Without long ranged dipole-dipole interac¬ 
tions, water far from the surface does not respond to the 
presence of this dipole layer, and £ (r) remains non-zero 
well into the bulk region. However, we can compensate 
for the averaged effects of the long ranged electrostatics 
through the introduction of the electrostatic LMF for an 
uncharged solute [13] 

V R (r) = J dr'p^(r')vi (|r - r'|), (6) 

where vi(r) = erf (r/a)/r is the long ranged, slowly vary¬ 
ing component of 1/r, separated with a smoothing length 
it = 4.5 A | herein, and in general a should be cho¬ 
sen to be greater than the nearest-neighbor distance in a 
fluid [13]. Inclusion of this renormalized solute potential 
in the GT water system leads to quantitative accuracy of 
both the electrostatic and orientational structure of in¬ 
terfacial water, evidenced by the curves labeled GT-LMF 
in Figure [2j 

In his seminal work on nonpolar solutes in aqueous so¬ 
lutions, Stillinger deduced that orienting an OH bond to¬ 
ward the interface provides the least energetic detriment 
to the hydrogen bond network of water [TJ . In GT water 
there are no opposing long ranged electrostatic interac¬ 
tions and the energetics of the hydrogen bond network 
alone determines the orientational preferences of water at 
the interface. However, this results in too high a proba¬ 
bility of pointing an OH bond toward the interface, illus¬ 
trating that while hydrogen bonding is a major driving 
force in determining the structure of water around large 
apolar solutes, it is not the sole determinant of the ob¬ 
served orientational preferences of interfacial water. 

In an earlier contribution, Stillinger and Ben-Naim 
initially postulated that the dipole and quadrupole mo¬ 
ments of water lead to a mean torque on a molecule at the 
interface with its vapor that orients the dipole moment of 
an interfacial water molecule toward the bulk liquid [23]. 
This behavior is reflected in the change of P (9^) upon the 
inclusion of long ranged interactions through V R , which 
provides the slowly-varying torque necessary to slightly 
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FIG. 2. (a) Polarization potential <3>(r) and (b) the corre¬ 
sponding electric fields £(r) obtained for a solute of I?hs ~ 
20 A in SPC/E and GT water, as well as for GT water in the 
presence of the electrostatic LMF (GT-LMF). (c) Probability 
distributions of the angle formed by the OH bond vector and 
the vector connecting the oxygen site with the center of the 
solute (0oh), for molecules within 1 A of the solute surface for 
the three systems shown in (a). The analogous distributions 
for the dipolar angle 9 M are shown in (d). 


turn the molecular dipoles of interfacial water in the di¬ 
rection of the bulk and obtain the desired orientational 
structure, evidenced by the distributions P ( 9 M ) shown 
in Figure [2}1. Therefore, the orientational structure of 
water at extended hydrophobic surfaces is a result of a 
delicate balance of the energetics of the hydrogen bond 
network and the mutipolar interactions arising from the 
electrical asymmetry of a water molecule, with the former 
dominating. 


IV. THE RESPONSE OF INTERFACIAL 
WATER TO UNBALANCED FORCES 

In this section, we examine the response of short 
ranged reference systems around solutes of varying sizes 
to the presence of very strong unbalanced forces like those 
seen in reality only for very large solutes. This provides 
a stringent test of the stability of the hydrogen bond net¬ 
work around small solutes even when subjected to strong 
perturbations. In order to accomplish this task, we scale 
the long ranged LJ LMF determined for a large solute of 
radius I?hs ~ 20 A by its radius, and then rescale the 
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field to the desired solute size, i?HS, 


4>m (r-, A, Rhs) = A^ri | -Rhs 


-Rhs 


(7) 


where (r) is the slowly-varying renormalized portion 
of the LMF shown in the inset of Figure[TJ). Here the no¬ 
tation 0^ (r; -Rhs) indicates that the field 0^ is a func¬ 
tion of r and that it was determined when a solute of 
radius i?ns is fixed at the origin. The fictitious, rescaled 
LMF is indicated by 0 ri, and the coupling parameter A 
is used to further adjust the magnitude of this field. In ef¬ 
fect we have taken the large unbalanced LJ force around 
a large solute, which Figure [f]) shows is strong enough 
to significantly perturb the large scale density profile of 
GTRC water when corrected with LMF theory, and ar¬ 
tificially applied it to a small scale system like that in 
Figure [T}i with an intact local hydrogen bond network. 
This provides insight into the very different response in¬ 
terfaces around small and large hydrophobic solutes have 
to repulsive forces over a wide range of magnitude as A 
is varied, including exceptionally large unbalanced forces 
seen in reality only near large hydrophobic solutes. 

In order to quantify the response of water to strong 
unbalanced forces, we focus on the A-dependence of the 
average number of water molecules in the solute solvation 
shell, (IV(A))^ ri , as well as the corresponding response 
function 


X(A) = - 






d\ 


( 8 ) 


where (• • • )^ ri indicates that the ensemble average is per¬ 
formed in the presence of the field 0 ri (r; A, I?hs) ■ The 

function (IV(A))^ ri is calculated for distances r < r min , 
where r m ; n is defined as the distance at which the den¬ 
sity distribution in the absence of the field reaches its 
first minimum. 

In the large scale hydration regime the broken hy¬ 
drogen bonds in the interfacial region effectively per¬ 
mit the interface to detatch from the solute and the in¬ 
terface is “soft” and fluctuating. We expect water to 
have a response qualitatively similar to that of simple 
liquids where drying occurs with increasing strength of 
0 ri- However, in the small length scale limit, while 
network fluctuations certainly occur, the hydrogen bond 
network is basically maintained around the solute. We 
thus expect that the small scale solute-water interface is 
“stiff” and highly resistant to perturbations unless they 
are strong enough to break hydrogen bonds. This should 
lead to behavior that is fundamentally different from that 
of a simple LJ fluid, which lacks such strong, local inter¬ 
actions. 

As postulated above, in the large solute regime, the 
behavior of (IV(A))^ ri and y(A) are qualitatively similar 
for both GTRC water and the WCA fluid (Figure^ and 
Figure [3ji). Gradual dewetting is observed with increas¬ 
ing field strength, until no molecules are present in the 



FIG. 3. Average number of truncated water and LJ molecules 
in the first solvation shell as a function of the coupling param¬ 
eter A for solutes of radii .Rhs ~ 3 A (a) and .Rhs ~ 15 A (b). 
Results are shown for both GTRC water and the WCA fluid 
(with the LJ parameters of SPC/E water), and (IV(A))^; has 
been normalized by its value in the case of zero field in order 
to make comparisons between the two fluids. The correspond¬ 
ing response functions are shown in (c) and (d), respectively. 
Solid lines in (a) and (b) are spline fits to (IV(A))^ and those 
in (c) and (d) are the negative derivatives of the corresponding 
fits. 


solvation shell region at high values of the coupling pa¬ 
rameter. In fact, as A is increased, a peak in the response 
function x is observed, indicative of a drying transition in 
the hydration shell of the solute; the details of the transi¬ 
tion differ between GTRC water and the WCA fluid due 
to differences in state points and interaction potentials. 

In the small solute regime, the WCA fluid displays 
signatures of a drying transition completely analogous to 
those seen in the large solute case with a simple shift in A. 
GTRC water, on the other hand, does not display char¬ 
acteristics of such nanoscale dewetting (Figure §■ and 
Figure[3b); (_/V(A))^ ri stays roughly constant and the re¬ 
sponse iunction fluctuates about zero. Using a typical 
geometric definition of a hydrogen bond EH22, we find 
that the average number of hydrogen bonds per molecule, 
for waters located between the solute and the position of 
the first maximum in the corresponding p(r), fluctuates 
around 3.5 for all A > 0, very close to the bulk value 
of 3.6 hydrogen bonds per water molecule. Therefore, 
the hydrogen bond network is maintained around the 
small solute for all studied values of A, and the strong 
local interactions of the hydrogen bond network prohibit 
drying at the solute surface, even in the presence of the 
extremely large external fields considered herein. 

The above-described results indicate that the underly¬ 
ing physics behind the solvation behavior in a LJ fluid is 
qualitatively similar in the small and large length scale 
regimes, dependent only on the magnitude of the unbal- 
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ancing potential arising from the bulk, while that of wa¬ 
ter qualitatively differs in the two regimes. In the large 
length scale regime, water behaves in a manner similar 
to a LJ fluid, with the unbalanced LJ attractions hav¬ 
ing a substantial impact on the solvation structure. For 
solutes smaller than the crossover radius, however, wa¬ 
ter wets the surface of the solute even in the presence 
of extremely large (though fictitious) unbalancing poten¬ 
tials; the hydration shell remains intact due to the great 
strength of the local hydrogen bond network. Therefore, 
interfacial fluctuations and the physics dictating where 
the length scale transition occurs is different for water 
than for simple, non-hydrogen bonding fluids. 


V. HYDROGEN BONDING SETS THE SCALE 
FOR THE CROSSOVER IN HYDRATION 
THERMODYNAMICS 

The above-described physical balance between hydro¬ 
gen bonding and interfacial unbalancing potentials also 
plays a key role in the solvation thermodynamics of ap- 
olar solutes. Gibbs free energies of solvation, AG, were 
calculated by performing equilibrium simulations of so¬ 
lutes with effective hard sphere radii Rhs < 13 A in incre¬ 
ments of AI?hs ~ 0.5 A. Due to poor phase space over¬ 
lap between neighboring windows, AAhs was decreased 
to 0.25 A to determine AG for solutes with i?ns > 7 A 
solvated by GTRC water. The solvation free energies 
presented herein were calculated using the Bennett ac¬ 
ceptance ratio or BAR [25] [26] method. To emphasize 
the crossover in the scaling behavior of the solvation free 
energies, we normalize AG by the surface area of the 
apolar solute (Figure [ 4 ]), AG = AG/47rAp S . 

In the small solute regime, Aps < Rc ~ 5.0 A, the hy¬ 
dration free energies are in agreement for all three mod¬ 
els. This illustrates that the hydration thermodynamics 
of small, nonpolar solutes are dictated by the local struc¬ 
ture of water alone, as would be expected from the con¬ 
clusions drawn above regarding solvation in the SPC/E, 
GT, and GTRC models. Indeed, the dominant role of lo¬ 
cal structure in the small solute regime is not restricted 
to water, as indicated by the agreement of the solvation 
free energies for LJ and WCA fluids for small solute sizes 
shown in Figure |4 Jd. 

The free energy for large solutes scales with surface 
area in both SPC/E water and the LJ fluid, and here long 
ranged interactions become increasingly important. Only 
small differences in AG are observed between SPC/E 
and GT water, reflecting the relatively small role of long 
ranged electrostatics in hydrophobic hydration [9|. LJ 
attractions, on the other hand, make a substantial con¬ 
tribution to the hydration free energy. Indeed because of 
the absence of these attractions, GTRC water completely 
lacks the plateau in AG for large solute sizes. 

The behavior of the GTRC water model can be ex¬ 
plained by noting that in the large solute regime, AG ~ 
PVs+jAs, where Vs and A$ are the volume and surface 



FIG. 4. Solvation free energies of apolar spheres per unit so¬ 
lute area as a function of solute radius, scaled by the effective 
diameter of a corresponding solvent molecule (d e ff = 2.75 A 
for water), for (a) SPC/E, GT, and GTRC water models, as 
well as (b) a LJ fluid and its corresponding WCA reference 
system. Error bars are smaller than the symbols shown. 


area of the solute, respectively, P is the pressure of the 
system, and 7 is the solute-water surface tension. In or¬ 
der to obtain the same bulk density as SPC/E water at a 
pressure of 1 atm, the GTRC model must be maintained 
at a pressure of roughly 3 katm. At this state the GTRC 
water model is far from liquid-vapor coexistence, and the 
pressure is large enough to make the PVs term dominate 
the behavior of AG for large solutes. 

However we have previously shown that GTRC water 
can indeed have a self-maintained liquid-vapor interface, 
but at a lower bulk density close to that of ice. The inter¬ 
face is maintained by the strong short ranged Coulomb 
attractive forces between donor and acceptor sites and 
the need to preserve as many hydrogen bonds as possi¬ 
ble [9]. However, because there are no unbalanced forces 
from LJ attractions, the surface tension is much smaller 
than that of the full SPC/E model. 

As shown in the curve labeled ‘GTRC-coex’ in Fig¬ 
ure [4^,, the solvation free energies in GTRC water near 
coexistence in both the small and large solute regimes 
are smaller in magnitude than those in SPC/E water. 
However it exhibits essentially the same crossover radius 
as the full SPC/E model and scales with solute surface 
area for large solutes. The behavior of AG below the 
crossover radius can be understood from our previous 
results for the bulk structure of the GTRC model near 
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coexistence |2j- The bulk coexistence density is close to 
that of ice and the hydrogen bond network has a more 
ordered tetrahedral structure that can more readily ac¬ 
commodate the formation of a cavity than is the case for 
SPC/E water. 

Although the solvation free energies of apolar solutes 
in water and in the LJ fluid exhibit qualitatively similar 
crossover behavior, they differ in one important respect: 
the length scale at which the crossover in solvation be¬ 
havior occurs. For the LJ fluid, the crossover radius is 
approximately equal to the diameter of a solvent parti¬ 
cle. At this solute size, the unbalanced forces from the 
LJ attractions of the bulk region become large enough to 
“pull” particles away from the solute surface, leading to 
drying. 

Although unbalanced LJ forces also exist when apolar 
particles of similar size are solvated by water, the possible 
disruption of strong local hydrogen bonds between inter¬ 
facial water molecules dominates the energetics, and the 
crossover occurs only when water is not able to maintain 
this network. This leads to an estimate for the crossover 
radius, Rc ~ 5 A, almost twice the diameter of a water 
molecule (2.75 A) and significantly larger than that found 
in a LJ fluid. As shown above, hydrophobic solvation in 
GTRC water near coexistence also displays a crossover 
in its scaling behavior at a value of Rc essentially the 
same as that of the full SPC/E model. Because GTRC 
water accounts only for the hydrogen bond network, we 
can conclusively say that the crossover in solvation be¬ 
havior is determined by the hydrogen bond network of 
water alone, occurring when the solute size is increased 
to a point beyond which it is impossible for this network 
to remain intact, consistent with the original arguments 
of Stillinger [33. 

Given the importance of the hydrogen bond network 
for small scale solvation in water, how can we rational¬ 
ize the success of the LCW theory [5 and related lat¬ 
tice models incorporating similar physics [35l 136] , which 
lack an explicit description of hydrogen bonds? These 
theories correctly describe the small scale physics driven 
by Gaussian density fluctuations in the bulk solvent and 
the large scale physics dominated by the formation of a 
vapor-like interface around a large repulsive solute. Ef¬ 
fective parameters controlling the transition between the 
two regimes are fit to experimental data for each partic¬ 
ular solvent. 

The key experimental parameters determining the 
transition length scale in the LCW theory are the liquid- 
vapor surface tension, and the bulk density and com¬ 
pressibility. The small compressibility and large surface 
tension of water compared to a LJ fluid implicitly ac¬ 
counts for the strength of the hydrogen bond network in 
bulk water and the difficulty of disrupting it by interface 
formation for large solutes. This allows the LCW theory 
to qualitatively describe the different transition length 
scales in both water and a LJ fluid [7j using the same 
basic framework. But LCW theory uses mean field ideas 
and square gradient and other approximations, and er¬ 


rors are seen in its detailed predictions for certain other 
properties like the interface width J7] . More detailed ap¬ 
proaches describing structure and fluctuations in both 
small and large length scale regimes are needed for quan¬ 
titative calculations. 

More recent work by Rajamani, Truskett, and 
Garde [3B] has clarified the relation between bulk ther¬ 
modynamics and the crossover radius. They suggested 
that the crossover radius is proportional to the Egelstaff- 
Widom length scale Zew = 7 «r, the product of the 
liquid-vapor surface tension 7 and the isothermal com¬ 
pressibility Kt m- Quantitative agreement can be 
achieved by using a microscopic compressibility that 
depends on the solute volume rather than the long 
wavelength bulk compressibility in conjunction with the 
solute-water interfacial tension to estimate the crossover 
radius Rc- 

A simple but stringent test of this idea is to compare 
the Egelstaff-Widom length scale of GTRC water near 
liquid-vapor coexistence to that of SPC/E water. As 
discussed above, the crossover radius in GTRC water is 
essentially the same as in SPC/E water. This is eas¬ 
ily rationalized from our microscopic understanding of 
the very similar behavior of the hydrogen bond network 
around the solute in GTRC and SPC/E water. If this 
simple physics is reflected in the Egelstaff-Widom length 
scale, this too should be nearly the same although both 
the surface tension and bulk compressibility differ con¬ 
siderably in the two models. 

Indeed, the compressibility k gtrc of GTRC water at 
T = 300 K and a pressure of 1 atm is 0.087 katm -1 , 
roughly a factor of two larger than that of SPC/E at the 
same state point, 0.045 katm -1 , while the surface ten¬ 
sion of the GTRC model 7 GTRC ss 27 mN/m, is about 
half of that of the SPC/E model 7 SPC / E ss 54.7 mN/m. 
Here the value for SPC/E water was taken from the work 
of Sedlmeier and Netz [39] and the surface tension of 
GTRC water was estimated by extrapolating the solva¬ 
tion free energies AG(i?Hs) presented in Section V to 
the limit Rhs —► 00. Thus, the Egelstaff-Widom length 
scales of SPC/E and GTRC water are nearly equal, 
C C/E = 0-24 A and Z ^™ 0 = 0-23 A, respectively, as 
expected. 


VI. ENTROPY CONVERGENCE IS A 
CONSEQUENCE OF THE HYDROGEN BOND 
NETWORK 

The temperature dependence of hydrophobic hydra¬ 
tion also displays features distinct from solvation in typ¬ 
ical van der Waals liquids. Specifically, hydration free 
energies A G of small apolar particles increase with in¬ 
creasing temperature along a significant portion of the 
coexistence curve until a maximum is reached. Above 
this temperature, free energies of solvation decrease with 
increasing temperature, a behavior typical of most flu¬ 
ids. Associated with this region of anomalous solvation 
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FIG. 5. Hard sphere solvation free energy AG per unit solute 
volume Vs as a function of temperature in (a) SPC/E and 
(b) GTRC water. The corresponding entropies of solvation 
AS as a function of T are shown in (c) and (d), respectively. 
Hard sphere radii are indicated in the legend. Solvation free 
energies as a function of solute size for T = 300 K, 325 K, 
350 K, 375 K, and 400 K are shown in the insets of (c) and (d) 
for the SPC/E and GTRC models, respectively. The arrows 
point in the direction of increasing temperature. 


is the phenomenon of entropy convergence , in which the 
hydration entropies, AS = — (<9A G/dT)p, intersect near 
a temperature of 400 K for a large range of solute sizes, 
although the location of the hydration free energy max¬ 
imum varies somewhat with solute size. Analogous to 
the discussion of the crossover length scale, the expla¬ 
nation of entropy convergence typically uses thermody¬ 
namic arguments, citing the small and nearly constant 
compressibility of water along the liquid-vapor coexis¬ 
tence line, relative to organic solvents piTil - l l 2 j . although 
explanations exist that do not hinge on the relative in¬ 
compressibility of bulk water [33! • 

In this section, we show that entropy convergence in 
water arises from the hydrogen bond network through its 
impact on bulk thermodynamics by studying the temper¬ 
ature dependence of hard sphere solvation in the SPC/E 
and GTRC water models near liquid-vapor coexistence. 
Simulations of bulk SPC/E and GTRC water were car¬ 
ried out at a pressure of 1 atm and temperatures rang¬ 
ing from 275K-500K and 225K-425K, respectively. Hard 
sphere solvation free energies in the small solute regime 
were determined by assuming Gaussian bulk density fluc¬ 
tuations mmm 


k B Tpl(T)V | 
2vvs( T ) 


k B T 

2 


In [ 27 rcry s (T)] , 


(9) 


where cry s = ((dTV) 2 )^ is the mean squared fluctuation 
in the number of molecules TV in a solute-sized probe 


volume V 3 , with SN = N — (N) Vs , and we consider the 
volume V s = 47ri?y S /3 of a spherical solute of radius 
Rhs herein. These solvation free energies were then fit 
to A G(T) = a + bT — cT 2 , and are plotted as lines in 
Figures [5ji and JHJr. Solvation entropies were determined 
from the negative derivative of these fits, and are shown 
in Figures [5jc and[5}l. 

The temperature dependence of hard sphere solvation 
is qualitatively similar in both SPC/E and GTRC wa¬ 
ter. In fact, entropy convergence is observed in the 
GTRC model, albeit at a convergence temperature T 
approximately 100 K less than the convergence tem¬ 
perature in SPC/E water; Xspc/e = 387 ±8 K and 
Pgtrc = 291 ±7 K, obtained from linear fitting of AS as 
a function of the heat capacity of solvation, A Cp{T) = 
T(dAS/dT)p, for several temperatures [33]. Despite this 
quantitative distinction, the fact that the minimal refer¬ 
ence network of GTRC water captures the phenomena 
of entropy convergence explicitly demonstrates that this 
signature of hydrophobic hydration is directly related to 
the energetics of the hydrogen bond network over a wide 
range of temperatures. 

Previous work has shown that the logarithmic term in 
Equation [9] has merely a secondary effect on entropy con¬ 
vergence, shifting T to somewhat lower values and A S(T) 
from zero to negative values [30]. Therefore, in order to 
obtain a qualitative, microscopic explanation for entropy 
convergence, we can neglect this term in the Gaussian ap¬ 
proximation for the free energy, and write the solvation 
entropy as 


A 

V 2 crys 


pl(T) [1-2 Ta P (T)} 


( 10 ) 


where ctp = — (d\np B /dT)p is the thermal expansion 
coefficient at constant pressure, which was determined 
by fitting the bulk densities to Laurent polynomials [45] . 
Here we have also assumed that the temperature depen¬ 
dence of the variance ay s can be neglected, as has been 
previously established m ■ Thus within the accuracy of 
Equation 10 entropy convergence is seen for A S(T) = 0, 


and an estimate of the convergence temperature can be 
obtained from the intersection of ap(T) and (2T) _1 . The 
convergence temperatures obtained for the SPC/E and 
GTRC models from Equation [lO] are roughly 420 K and 
330 K, respectively, in reasonably good agreement with 
the results presented above, although T will always be 
overestimated in this approximation. Nonetheless, the 
difference between the convergence temperatures of the 
two models is quantitatively captured by this estimation, 
indicating that additional T-dependences arising in AG 
are similar in the two models, and these have been dis¬ 
cussed in detail elsewhere [40] [42] . 

In this simplified Gaussian framework the behavior of 
ap(T) plays a key role in entropy convergence. In the 
case of SPC/E water, the thermal expansion coefficient 
vanishes at the temperature of maximum density near 
248 K P [45]. As shown in Figure [6j 2 Tap(T) then 
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FIG. 6. Thermal expansion coefficient multiplied by twice the 
temperature (left axis, closed symbols) and average number of 
hydrogen bonds per water molecule (right axis, open symbols) 
for SPC/E and GTRC water. 







R( A) 


FIG. 7. Potential of mean force, W(r), between two UA 
methane particles in SPC/E, GT, and GTRC water. 


increases with increasing temperature but remains less 
than one until about 420 K, where entropy convergence 
is predicted to occur. The thermal expansion coefficient 
of GTRC water behaves in a qualitatively similar man¬ 
ner with 2 Tap(T) remaining less than one until about 
330K, although ap(T) is never negative, because this 
model lacks a density maximum near liquid-vapor coex¬ 
istence [9j. The behavior of the thermal expansion coef¬ 
ficient is a direct consequence of the energetics of the 13- 
bond network in both models. At ambient temperatures, 
the average number of hydrogen bonds per molecule ap¬ 
proaches four in both SPC/E and GTRC water [9]. With 
increasing temperature, thermal fluctuations increasingly 
disrupt the entropically unfavorable hydrogen bond net¬ 
work in both models (Figure |b|, which leads to an in¬ 
crease in the thermal expansion coefficient. However, the 
lower density of GTRC water permits more fluctuations 
as the temperature is increased, consistent with its larger 
compressibility and a more rapid increase in ap(T), lead¬ 
ing to a lower convergence temperature. 

We also determined the temperature dependence of 
large solute solvation free energies following the descrip¬ 
tion in the previous section. After the length scale tran¬ 
sition, solvation is dominated by interfacial physics. As 
evidenced by the insets in Figure [5] hard sphere solvation 
free energies in this regime decrease with increasing tem¬ 
perature for both models, following the T-dependence of 
the surface tension, just as is the case for LJ solvation. 


VII. LONG RANGED INTERACTIONS AND 
THE SIZE DEPENDENCE OF HYDROPHOBIC 
ASSOCIATION 

In this section, we examine the role of the various short 
and long ranged forces in the thermodynamics of hy¬ 
drophobic association. In order to accomplish this task, 


we consider the association of pairs of spherical solutes, 
one pair in which both solutes are in the small-scale 
regime, while the other pair consists of two large solutes. 
We first examine the free energy as a function of solute- 
solute distance, R , 

/3W(R) = -lnP(R), (11) 

where P(R) was obtained by umbrella sampling with the 
harmonic biasing potential 

U h UR) = ^(R-R*) 2 , (12) 

R* is the desired value of R. and k is a force con¬ 
stant tuned to achieve adequate overlap between neigh¬ 
boring windows. The probability distribution P(R) 
was then constructed from the set of biased simulations 
using the multistate Bennet acceptance ratio method 
(MBAR) [27]. 

We first focus on hydrophobic association in the small 
scale regime, and consider the association of two united 
atom (UA) methane models, which are simply LJ par¬ 
ticles with length and energy parameters of dMe-Me = 
3.73 A and CMe-Me = 1-234 kJ/mol, respectively 281 . 
Methane-water interactions were obtained from Lorentz- 
Berthelot mixing rules. 

The potentials of mean force, W(R), shown in Figure[7] 
for the association of two UA methanes are nearly identi¬ 
cal for all water models under consideration. Therefore, 
not only does the hydrogen bond network dictate the 
solvation structure around individual small solutes, but 
also the association of solutes in this length scale regime, 
as expected from the results presented in the previous 
sections. 

We now consider the association of two large C60 
fullerene molecules in the various models of water. Each 
C60 is represented as a single site using the coarse- 
graining procedure prescribed by Girifalco Earn], such 
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that the fullerene-fullerene interaction is given by 


U FF (R) = —a 

+ c 


1 


1 


s(s - l) 3 
1 


s(s +1) 3 
1 


s(s — l) 9 s(s + l) 9 s 10 


(13) 


where a = 4.4775 kJ/mol, £ = 0.0081 kJ/mol, s = R/2rj, 
and r/ = 3.55 A. The C60-water interaction potential is 


U wF (r) 


4Ve wF ^ 

rrj 


0wF 
rj — r 


10 



1 

8 


( Q~wF V 

\V - r) 


f 0\vF \ 

\V + r J 


4 ' 


(14) 


where N = 60, cr wF = 3.19 A, and e wF = 0.392 kJ/mol. 
Previous work has shown that this coarse-grained water- 
C60 interaction provides a very good representation of 
the solvation structure in the corresponding atomically- 
detailed water-C60 system m- 

The water-C60 interaction potential U wF (r) leads to 
a hydrophilic particle due to the high density of carbon 
sites on the surface of the C60 molecule. Therefore, we 
also consider a hydrophobic particle obtained by using 
only the repulsive water-C60 and C60-C60 forces. This is 
obtained by performing a WCA-like separation of the po¬ 
tentials U FF and U wF to obtain the corresponding purely 
repulsive potentials [7o,ff and t/o, W F, as detailed above 
for f/sw ■ 

We further separate the potential of mean force as 
W(R) = W vac (R) + W solv (R), where W vac (R) and 
Wsoiv(-R) are the vacuum and solvent-induced portions 
of the PMF, respectively, focusing on the latter contribu¬ 
tion herein. The solvent-induced PMFs between purely 
repulsive C60 particles in the SPC/E and GT water mod¬ 
els, shown in Figure[8ji, are indicative of the hydrophobic 
effect; the association of two large apolar particles in wa¬ 
ter is barrierless, although the free energy of association 
is slightly lower in GT water due to its lower surface ten¬ 
sion. Previous work has shown that the collapse of two 
extended hydrophobic surfaces proceeds by the formation 
of a vapor tube EH 132, in which solvent molecules are 
evacuated from a cylindrical region between the two hy¬ 
drophobes, and we will show below that the association 
of two repulsive fullerenes also occurs by this mechanism. 

In GTRC water, however, the PMF kF so i v (-R) displays 
a slight barrier at R ss 15 A, and another significantly 
higher barrier at R ~ 12 A, as shown in Figure [8^i. Be¬ 
cause the C60-C60 distance does not explicitly account 
for changes in the behavior of the aqueous solvent, it is 
not a good reaction coordinate to study the association 
of two large hydrophobes on its own .31, 32] and W (R) 
cannot provide an explanation for the appearance of this 
barrier in W so i v (i?). 

To understand hydrophobic association in GT and 
GTRC water, we calculate the free energy as a function of 
the C60-C60 distance R and the density p v of water in a 
cylindrical volume of radius 3.75 A between the particles. 




FIG. 8. Solvent-induced potential of mean force, II4oiv(A), 
between two (a) purely repulsive and (b) attractive coarse¬ 
grained C60 particles in SPC/E, GT, and GTRC water. In¬ 
sets in (a) and (b) show the corresponding nonuniform den¬ 
sities around a single coarse-grained C60 immersed in each 
water model. The gray curve in (b) is the GTRC PMF from 
panel (a). 


This two-dimensional free energy landscape is given by 
/ 3W(R,p v ) = —In P(R,p v ), where P(R,p v ) was calcu¬ 
lated using the indirect umbrella sampling method [33] 
to bias the number of particles in the volume v. The 
harmonic potential in Eq. |T2] was used to bias R. Again 
MBAR was used to reconstruct the probability distribu¬ 
tion from these biased simulations [27] . 

The free energy surface shown in the top panel of Fig¬ 
ure [9] indicates that hydrophobic collapse in GT water 
(or SPC/E water) is indeed driven by the barrier-less 
formation of a vapor tube (31j [32] at a C60-C60 distance 
between 14 and 15 A. Hydrophobic collapse in GTRC 
water, on the other hand, does not follow this mecha¬ 
nism because capillary evaporation in the inter-fullerene 
region has been suppressed by the removal of LJ attrac¬ 
tions in the solvent. This is consistent with the lack of 
drying at the interface of a single repulsive solute, as ev¬ 
idenced by the nonuniform densities shown in the inset 
of Figure [Hjj, and would be anticipated from the results 
presented in Section IV. 

Instead, the free energy minimum in GTRC water 
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FIG. 9. Free energy as a function of C60-C60 distance, R, and 
density of water in the observation volume v with respect to 
that in the bulk, p v /pB, for the association of two hydrophobic 
fullerene particles in (a) GT and (b) GTRC water models. 
Contour lines are spaced in increments of k^T. 


(for a specific value of R ) remains at liquid-like densi¬ 
ties as the C60-C60 distance is decreased, until the water 
molecules cannot physically remain between the fullerene 
particles due to repulsive core overlap near R ss 12 A. 
Only at this point are the solvation shell water molecules 
in the inter-fullerene region expelled. This expulsion of 
water molecules in the observation volume causes the 
large free energy barrier observed at the same inter- 
fullerene distance in the one-dimensional Wsoiv(-R) for 
GTRC water shown in Figure [8ji. 

Instead of artificially suppressing capillary evaporation 
between large hydrophobes by removal of solvent LJ at¬ 
tractions as in GTRC water, we can directly counteract 
the unbalanced LJ interfacial forces leading to evapora¬ 
tion in the GT or full water models by making the solutes 
sufficiently hydrophilic. LMF theory would predict very 
similar behavior for these two systems. This is accom¬ 
plished by using the full {/ff and U w f potentials to de¬ 


scribe fullerene-fullerene and water-fullerene interactions, 
respectively. Inclusion of the water-C60 attractive inter¬ 
actions leads to an almost perfect cancellation of these 
unbalanced forces, as evidenced by the good agreement 
of the SPC/E and GT nonuniform densities with that of 
the GTRC model, shown in the inset of Figure [HJo. 

These strong solute-water attractions, arising from the 
high surface density of carbon atoms, render the C60 
molecule hydrophilic, and the associated solvent-induced 
PMFs are repulsive for all distances. This indicates that 
water opposes the association of two such particles, in 
accord with previous results [33]. Because of the effec¬ 
tive hydrophilicity of the particles, capillary evaporation 
between the particles does not occur, and W so iv(.R) is the 
same for all three models for R > 12 A. At smaller sepa¬ 
rations water is forcibly expelled from the inter-fullerene 
region due to overlap with the repulsive cores of the so¬ 
lutes and then differences arise due to the differing pres¬ 
sure of the systems. 

The two-dimensional PMF W{R,p v ) was also calcu¬ 
lated for the case of hydrophilic fullerene particles in GT 
water, and is shown in Figure [lO] This PMF is qualita¬ 
tively very similar to that shown for hydrophobic collapse 
in GTRC water in Figure [9 ]d as expected. As R is de¬ 
creased, the free energy minimum as a function of p v 
remains in regions of liquid-like densities. It is not until 
very small R, less than 12 A, that W{R 1 p v ) develops a 
minimum at low p v , indicating a global free energy min¬ 
imum at the contact state. In fact, the solvent induced 
PMF Wsoiv(i?) between hydrophobic solutes in GTRC 
water is nearly identical to the PMFs obtained between 
hydrophilic solutes in all models until water is expelled 
from the inter-fullerene region, R < 12 A, as illustrated 
by the curve labeled l GTRC-Hphob’ in Figure [8 Jd. In 
contrast to what is found for the association of large hy¬ 
drophobic particles, the solvent opposes association and 
the contact state is stabilized by the large solute-solute 
attractions between hydrophilic fullerenes. 


VIII. CONCLUSIONS 

We have used short ranged variants of the SPC/E wa¬ 
ter model [5J in conjunction with LMF theory to examine 
the crossover in the behavior of hydrophobic hydration 
with increasing solute size. While small scale solvation 
is determined exclusively by the local structure of wa¬ 
ter, i.e. the hydrogen bond network, long ranged inter¬ 
actions are important for the accurate description of the 
hydration of large apolar solutes. Dispersion interactions 
lead to the phenomena of drying at extended hydropho¬ 
bic interfaces, while long ranged dipolar interactions are 
essential for the description of the orientational ordering 
of water in the vicinity of a large solute, as well as for 
interfacial electrostatic properties. 

The truncated GT and GTRC water models also pro¬ 
vide insight into hydrophobic interactions between so¬ 
lutes in the small and large length scale regimes. The 
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FIG. 10. Free energy as a function of C60-C60 distance, R, 
and density density of water in the observation volume v with 
respect to that in the bulk, p v /pB, for the association of two 
hydrophilic fullerene particles in GT water. Contour lines are 
spaced in increments of k^T. 

local structure of water, dictated by the hydrogen bond 
network, is found to govern the association of two small 
scale solutes, a concept which has been successfully ex¬ 
ploited to provide a theoretical framework for describing 
hydrophobic hydration and association at small length 
scales [33] . Moreover, previous work has shown that 
coarse-grained models, whereby water molecules inter¬ 
act via a single spherically symmetric pairwise potential, 


can reproduce the thermodynamics of association of two 
methanes [461 [47]. From the results presented here, it is 
not surprising that such coarse-grained models can cap¬ 
ture features of small scale hydrophobicity, since these 
models also describe the bulk structure of water with 
near quantitative accuracy. 

The association of two large scale hydrophobes involves 
the formation of an inter-solute vapor tube, and the un¬ 
balanced forces arising from water-water LJ attractions 
are found to be of the utmost importance for this mech¬ 
anism of hydrophobic association. In this regime the 
coarse-grained water models will fail completely. Cancel¬ 
lation of the effects of interfacial unbalanced forces, ei¬ 
ther by explicit removal of solvent-solvent LJ attractions 
(as in GTRC water) or by addition of large solute-water 
attractions that counterbalance these forces, suppresses 
capillary evaporation between two large solutes. As a 
result the solute surface is wet by the aqueous solvent, 
and free energy barriers to the association of two large 
hydrophilic solutes exist. In all these cases comparison 
of results in the full model with those from the short- 
ranged GT and GTRC water models provides a simple 
and physically suggestive way to disentangle the effects of 
longer ranged dispersive and Coulomb interactions from 
properties of the local hydrogen bond network. 
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